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^ ■ Abstract 

O 
O 

Q_i| It has recently been pointed out that viable sparticle mass spectra can 

q_|' be generated in Yukawa unified 5*0(10) super symmetric grand unified models 

consistent with radiative breaking of electroweak symmetry. Model solutions 
are obtained only if tan (3 ~ 50, [i < and positive D-term contributions to 
scalar masses from SO (10) gauge symmetry breaking are used. In this paper, 
we attempt to systematize the parameter space regions where solutions are 
obtained. We go on to calculate the relic density of neutralinos as a function 
of parameter space. No regions of the parameter space explored were actually 
cosmologically excluded, and very reasonable relic densities were found in 
much of parameter space. Direct neutralino detection rates could exceed 1 
event /kg/day for a 73 Ge detector, for low values of GUT scale gaugino mass 
m i/2- We also calculate the branching fraction for b —* sj decays, and find 
that it is beyond the 95% CL experimental limits in much, but not all, of the 
parameter space regions explored. However, recent claims have been made 
that NLO effects can reverse the signs of certain amplitudes in the b — > S7 
calculation, leading to agreement between theory and experiment in Yukawa 
unified SUSY models. For the Fermilab Tevatron collider, significant regions 
of parameter space can be explored via bbA and bbH searches. There also exist 
some limited regions of parameter space where a trilepton signal can be seen 
at TeV33. Finally, there exist significant regions of parameter space where 
direct detection of bottom squark pair production can be made, especially for 
large negative values of the GUT parameter Aq. 
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I. INTRODUCTION 



Supersymmetric SO (10) grand unified theories |j] (GUTS) are highly motivated choices 
for particle physics models which encompass and go beyond the Standard Model (SM). In 
particular: 

• All three forces in the SM are unified by a simple Lie group. Since the model is 
supersymmetric, the hierarchy between weak and GUT scales is naturally stabilized. 

• The fifteen distinct fermions of each SM generation plus a SM gauge singlet right- 
handed neutrino state are economically included in the sixteen dimensional spinorial 
representation of 5*0(10). 

• The group 5*0(10) is anomaly-free, and provides an explanation for cancellation of 
triangle anomalies that are otherwise rather ad-hoc in both the MSSM and SU (5) 
SUSY-GUT models. 

• The SM weak hypercharge assignments can be derived from broken 50(10) using only 
the lowest dimensional 50(10) Higgs multiplets. This has been interpreted as group 
theoretic evidence of 50(10) grand unification [0. 

• The singlet neutrino superfield(s) N c may develop a Majorana mass at very high energy 
scales ~ 10 11 — 10 16 GeV. Combining these with the usual Dirac neutrino masses that 
can now develop, one is lead to (dominantly) left-handed neutrino masses at the eV 
scale (or below), while right handed neutrinos are pushed beyond observability, via 
the see-saw mechanism [[|. The resulting neutrino masses can easily be in accord 
with recent data from solar and atmospheric neutrino detection experiments, though 
a detailed model is necessary for an explanation of the observed mixing pattern. 

• In minimal 50(10), there exists only a single Yukawa coupling per generation, so that 
Yukawa couplings also unify at the GUT scale. This is especially constraining for 
the third generation. In addition, predictive 50(10) models can be created which 
encompass the fermion masses of the lighter generations as well ||. 

• In 50(10) models, i?-parity conservation is a natural consequence of the gauge sym- 
metry of the model ||. When 50(10) is broken, if the breaking occurs via certain 
"safe" representations of the Higgs fields, i?-parity is still conserved even in the weak 
scale theory ||. Along the same lines, it has been shown that in SUSY models with 
gauged B — L (this includes 50(10) models) symmetry and a renormalizable see-saw 
mechanism, i?-parity must be conserved exactly, even in the weak scale Lagrangian 

1- 

• In 50(10) models, baryogenesis in the early universe can be explained as a consequence 
of decays of (out of thermal equilibrium) intermediate scale right-handed neutrino 
states |5J. 

Motivated by these considerations, in this paper we investigate a variety of phenomeno- 
logical implications of Yukawa-unified 50(10) SUSY-GUT models. Throughout our study, 
we assume the following: 
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1. A minimal SO (10) SUSY GUT model with a single Yukawa coupling per generation, 
valid at energy scales Q > M GUT ~ 2 x 10 16 GeV. 

2. The gauge symmetry 50(10) breaks via some (to be specified) mechanism directly to 
the MSSM gauge group at Q = Mqut- We assume universality of soft SUSY breaking 
mass parameters at Q = Mqut- However, as a consequence of the reduction in rank 
of the gauge symmetry breaking, additional D-term mass contributions are generated 
for scalar masses. 

3. Electroweak symmetry is broken radiatively (REWSB). 

With these assumptions, the spectrum of supersymmetric particles and Higgs bosons can 
be calculated. We examine various phenomenological consequences of the sparticle mass 
spectrum generated. We do not address the related question of mass generation of first and 
second generation SM fermions; this issue has been investigated in a number of other papers 
Q. We also do not explore the specific mechanisms involved in gauge symmetry breaking; 
several studies along these lines are contained in Ref. M. 

A crucial issue is to find models with Yukawa coupling unification which reproduce the 
observed masses for the third generation fermions. One may assume Yukawa unification at 
Q = Mqut, and evolve gauge and Yukawa couplings and soft SUSY breaking (SSB) masses 
to the weak scale, and only accept solutions which generate third generation fermion masses 
in agreement with measurements within tolerances. Alternatively, one may begin with the 
central values of measured fermion masses at low energies, and accept only solutions which 
give Yukawa unification to within set tolerances at the GUT scale. In this paper, we adopt 
the latter bottom-up approach. A crucial aspect for this approach is to correctly calculate 
the weak scale third generation Yukawa couplings. Throughout our work, we fix pole masses 
m b = 4.9 GeV, m t = 175 GeV and m T = 1.78 GeV. 

The mass spectrum of SUSY particles in minimal supersymmetric 50(10) constrained 
by radiative electroweak symmetry breaking has been studied previously in a number of 



papers [ 10 — 20| . Unification of bottom, tau and top Yukawa couplings was found to occur 



at very large values of the parameter tan (3 ~ 50 — 60, and specific spectra were generated 



for values of m t ~ 190 GeV [11|. Assuming universality of soft SUSY breaking masses at 
Mqut, it was found [12], 14] that Yukawa unification consistent with radiative electroweak 
symmetry breaking could also occur for m t < 170 GeV as long as ^ 300 GeV. This 
generally leads to sparticle masses far beyond the reach of the CERN LEP2 or Fermilab 
Tevatron pp colliders. For values of m t ~ 175 GeV, solutions including radiative electroweak 
breaking were very difficult to achieve. In Ref. [TP]], the SUSY particle mass spectrum was 



investigated with non-universal SSB masses. Various solutions were found, but ad hoc non- 



universality conflicts with the 5O(10) symmetry. In Ref. ||20|1 , it was suggested that 5O(10) 
D-term contributions to scalar masses had the correct form to allow for successful radiative 
electroweak symmetry breaking and the computation of weak scale SUSY particle masses, 
but explicit solutions were not presented. 

In a previous paper [pl |, we have shown that Yukawa unified solutions to the super- 
particle mass spectrum could be generated (even for m t = 175 GeV) assuming universality 
within 5O(10), but only if additional D-term contributions to scalar masses were included. 
Solutions with Yukawa unification good to 5% and REWSB were found only for values of 
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/i < and positive values of the D-term contributions. A crucial ingredient in the calcula- 
tion is the inclusion of SUSY loop corrections to fermion masses [fLS , 12j, 22^1 , from which the 
third generation weak scale Yukawa couplings are calculated. For fi > 0, the diminution in 
the 6-quark Yukawa coupling is so severe that not even b — r unification can be acheived 
for any values of tan/? : 2 — 60. For values of tan/? > 60, the Yukawa couplings diverge 
during renormalization group running from the weak scale to the GUT scale, so that these 
values of tan/3 form an upper limit to the allowed parameter space. For fi < 0, assuming 
universality of SSB scalar masses, a much lower parameter space limit of tan (3 < 45 is found. 
For higher values of tan (3 we find fi 2 < 0, and REWSB falters essentially because the Higgs 
mass squared rn 2 Hd is driven to smaller weak scale values than m 2 Hu . However, if one allows 
positive D-term contributions to scalar masses, then m 2 Hu < rn 2 Hd already at Mqut, and the 
region of allowed tan f3 expands to tan (3 ~ 50 — 55. This is sufficient to allow for t — b — t 
Yukawa unification for many models with < 0. Of course, the D-term contributions must 
be added to the various squark and slepton soft masses as well, so that the D-terms leave a 
distinctive imprint on the superparticle mass spectrum 2^] . The most striking possibility is 
that the light bottom squark could be by far the lightest of all the squarks, and furthermore, 
even within range of Fermilab Tevatron searches. In addition, left- slepton masses could be 
lighter than right slepton masses, in contrast to expectations from the SUSY models with 
universality at Q — Mqut- This latter effect could certainly be extracted from precision 
measurements at linear e + e~ colliders operating at y/s > 1 TeV, and would be evidence in 
favor of a Yukawa unified SO (10) SUSY GUT model. 

In this paper, we address many of the observable aspects of Yukawa-unified SO (10) 



models in some detail. In Sec. II, we summarize the results from Ref. | 2l| , and try to 
systematize the parameter space expected in Yukawa unified SO (10) model framework. In 
Sec. Ill, we present calculations of the cosmological relic density of neutralinos in this class 
of models. Over much of the parameter space, we find very reasonable relic densities, unless 
■nig ~ ^j 4 , where s-channel annihilation through a Higgs resonance can lower Qh 2 to very 
tiny values. Given a reasonable relic density, it is well-known that rates for direct neutralino 
detection experiments are highest at large tan/3. We show direct detection rates in S'O(IO) 
parameter space, and show that much of the space should be explorable by neutralino- 
nucleon scattering experiments. In Sec. IV, we evaluate the rate for b — > 37 decays. Much 
of the parameter space is excluded by the measured value of the branching fraction for the 
decay b — > 57, but some distinctive regions remain viable. These results may be subject 
to modification if there is additional non-universality or inter-generational squark mixing 
at Q = Mgut-i or if CP violating phases are significant. In addition, recent claims have 
been made that NLO QCD corrections to b — > 57 loop diagrams can reverse the signs of 
certain amplitudes, leading to agreement between experimental measurements and theory 
predictions from Yukawa unified SUSY models. In Sec. V, we examine regions of parameter 
space accessible to searches for supersymmetry at the Fermilab Tevatron collider. We find 
three possibilities: i) signals from Abb or Hbb associated production, ii) W\Z 2 — > 3£ signals 
if tan f3 is not too big, and in) signals from direct production of b\bi pairs, where observable 
rates occur for large negative values of Aq. Experiments at the Large Hadron Collider (LHC) 
would, of course, decisively probe the model. We end in Sec. VI with a summary of our 
results, and some concluding remarks. 
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II. MINIMAL 50(10) MODEL PARAMETER SPACE 



In the minimal supersymmetric 50(10) GUT model, the fifteen matter superfields of 
each generation of the MSSM plus a gauge singlet neutrino superfield N c are included in 
the 16-dimensional spinorial representation ip of 50(10). In addition, the two Higgs doublet 
superfields H u and H d are embedded in a single 10-dimensional Higgs superfield (p. The 
superpotential includes the term 

/ 3 f#$4> + ■■■ 

responsible for quark and lepton masses, with / the single Yukawa coupling per generation 
in the GUT scale theory. The dots represent terms including for instance higher dimensional 
Higgs representations and interactions responsible for the breaking of 5O(10). We neglect 
Yukawa couplings and fermion masses for the first two generations. As usual, SUSY break- 
ing, is parametrized by associated soft SUSY breaking mass terms. We assume a common 
mass mis f° r all matter scalars and a mass mio for the Higgs scalars, along with a universal 
gaugino mass m^, and common trilinear and bilinear SSB masses A and B. Motivated by 
apparent gauge coupling unification in the MSSM, 5O(10) is assumed to break directly to 
the gauge group SU(3)c x SU{2) L x U{l)y at Q — M GUT , where the reduction in rank of 
the gauge symmetry leads to distinct D-term contributions to scalar field mass terms. The 
scalar field squared masses at Q = Mqut are then given by 

m Q= m E = m U = m 16 + M D 

m 2 D = m\ = m\ % — 3Mp 
m L, d = m ?oT2M|,, 



where Mp parametrizes the magnitude of the Z)-terms, and can, owing to our ignorance of 
the gauge symmetry breaking mechanism, be taken as a free parameter, with either positive 
or negative values. \Md\ is expected to be of order the weak scale. Thus, the model is 
characterized by the following free parameters 1 : 

77116, 772io, Mf,, 7771/2, A , sign(fl). 

The value of tan (3 will be restricted by the requirement of Yukawa coupling unification, and 
so is tightly constrained to a narrow range around tan/3 ~ 50. 

In our previous paper pi] , our procedure was as follows. We generated random samples 



of model parameter sets within the ranges, 



< mis < 1500 GeV, 
< mio < 1500 GeV, 



1 If the effective theory below the GUT scale is the MSSM plus a right-handed neutrino, then 
several other parameters enter the discussion. Typically, these have only a small effect on the 
SUSY particle mass spectrum: see Ref. [24| 
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< m 1/2 < 500 GeV, 
-500 2 GeV 2 < M 2 < +500 2 GeV 2 , (2.1) 

45 < tan (3 < 55, 
-3000 GeV < A < 3000 GeV and 
ji > or ji < 0. 

For each parameter set, we then calculated the non- universal scalar masses according to 
formulae given above, and enter the parameters into the computer program ISASUGRA. 



ISASUGRA is a part of the ISAJET package |25j which calculates an iterative solution to 
the 26 coupled RGEs of the MSSM. 

To calculate the values of the Yukawa couplings at scale Q = Mz, we began with the pole 
masses m& = 4.9 GeV and m T = 1.78 GeV. We calculated the corresponding running masses 
in the MS scheme, and evolve and m T up to Mz using 1-loop SM RGEs. At Q = Mz, we 
included the SUSY loop corrections to mj and m T using the approximate formulae of Pierce 
et al. |[22|| , and then obtained the corresponding Yukawa couplings. A similar procedure is 
used to calculate the top quark Yukawa coupling at scale Q = m t . We assume a pole mass 
m t = 175 GeV. 

Starting with the three gauge couplings and t, b and r Yukawa couplings of the MSSM at 
scale Q = M z (or m t ), ISASUGRA evolves the various couplings up in energy until the scale 
where g\ = g<i, which is identified as Mqjjt, is reached. The GUT scale boundary conditions 
are imposed, and the full set of 26 RGE's for gauge couplings, Yukawa couplings and relevant 
SSB masses are evolved down to Q ~ M weak , where the renormalization group improved 
one-loop effective potential is minimized at an optimized scale choice Q = y j™>t L ' m t R and 
radiative electroweak symmetry breaking is imposed. Using the new spectrum, the full set 
of SSB masses and couplings are evolved back up to Mqut including weak scale sparticle 
threshold corrections to gauge couplings. The process is repeated iteratively until a stable 
solution within tolerances is achieved. We accepted only solutions for which the Yukawa 
couplings f t , fb and f T unify to within 5%. This constraint effectively fixes the value of 
tan/? to a narrow range of 49 ± 3. Yukawa unified solutions are found only for values of 
H < 0. This latter result is in accord with findings of Pierce et al. [^] where solutions for 
b — r unification in models with universality were found only for one sign of \x when tan (3 
was large. We also require the lightest SUSY particle to be the lightest neutralino, and that 
electroweak symmetry is successfully broken radiatively. 

In this paper, we have used upgraded ISAJET 7.47 including 2-loop Yukawa coupling 
RGEs [BH] with weak scale threshold corrections P7f to check our previous solutions. The 



updated solutions as a function of model parameter space are presented in Fig. [1]; no 
apparent differences can be seen compared with our earlier results which used just 1-loop 
RGEs for Yukawa coupling evolution. The Yukawa unified solutions to the SUSY particle 
mass spectra are found to occur typically with m 16 < m 10 ^ 1.5m 16 and for 0.1m 16 < M D < 
.35rai6- Too small of a D-term will not yield enough splitting in the Higgs boson soft masses 
to yield REWSB, while too large a D-term will lead to third generation scalar field masses 
below experimental limits, or to charge and/or color breaking minima in the scalar potential. 
Our requirement of 5% unification is determined by defining the variables rv, and r tT , 
where for instance r br = max(fb/ f T , fr/ fb)- We then require R = max(r br , r tb , r tr ) < 1.05. 
The methodology employed in Ref. plj is useful for determining the general regions of 
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model parameter space where Yukawa-unified solutions with REWSB can be found. How- 
ever, in order to survey the predictions of model parameter space for observable quantities 
of experimental interest, it is useful to develop a more systematic approach to the model 
parameter space. For most of the results in this paper, we will take tan (3 = 50 as an indica- 
tive central value. We are also restricted to fi < 0. Guided by the results of Fig. |I]a, we can 
also adopt a central value mio = 1.25mi6, and from Fig. Mo = |mi 6 or Mr, = \mi§. 
Then, most of the variation in parameter space comes from the values of mi6 and mi/2, for 
which, from Fig. |l]c, there is no obvious correlation. 

We show in Fig. [2] the resulting m 16 vs. myz plane, taking as well A = 0, for a) 
the smaller and b) the larger value of Mj> The region shaded by solid dots is excluded 
on the left by negative values of m\, and on the right by negative values of /z 2 : in either 
case, the REWSB constraint is violated. The unshaded region gives viable solutions to 
the superparticle mass spectra, although Yukawa unification can vary throughout the plot 
from regions with R < 1.05 to R ~ 1.25. The value of R can be fine-tuned to less than 
1.05 throughout much of the parameter space plane by adjusting the value of tan/5 within 
the band in Fig. To gain an idea of the superparticle and Higgs boson mass spectrum 
throughout the parameter space plane, we show dashed contours of squark mass ma L = 1000 
(lower) and 2000 GeV (upper contour). As noted in Ref. |ZT|, the first and second generation 
squark masses are bounded below by about 700 GeV, and are usually much heavier in the 
Yukawa-unified 5*0(10) model. The chargino mass contours of = 150, 250 and 350 GeV 
are shown as solid contours, increasing with mu2- They asymptote towards the right-hand 
parameter space boundary, which is given by solutions where /x 2 = 0. For large values of the 
weak scale gaugino mass M2, |/^| ~ M2, so that the lightest chargino and neutralinos contain 
significant higgsino components. Although first and second generation squarks and sleptons 
are quite heavy and beyond the reach of LEP2 and Tevatron experiments, the charginos 
can be light, and within the reach of these collider facilities. Finally, the dotted contours 
show values of the pseudoscalar Higgs mass rriA = 150, 250 and 350 GeV, from left to right. 
Regions of parameter space with rriA — 100 — 200 GeV and large tan /3, as in this model, 
may have Higgs bosons accessible to Fermilab Tevatron collider searches via the bbA and 
bbH production mechanisms p8| . 

In Fig. 0, we show similar plots of model parameter space, this time taking Mr> = |mi6, 
but a) A = — mi6 and b) A = +mi 6 . The mass contours are the same as in Fig. ||. In 
Ref. |[21|| ) it was shown that the light bottom squark b\ could well be by far the lightest of 
all the squarks in Yukawa-unified SO (10). For negative A parameters, mixing effects in the 
bottom squark mass matrix can give rise to even smaller values of ~ 100 — 300 GeV, 
which may be accessible to direct searches at the Fermilab Tevatron P^. In addition, as 
can be seen, the range of parameter space is increased beyond the corresponding result for 
Aq = 0. Alternatively, taking Aq = +mi6 yields typically heavier bottom squark masses, 
and a smaller range of accessible parameter space. 

In Fig. £|, we show again the parameter space plane for M D = |mi 6 but with a) A Q = 
— m 16 and b) A = 0. The planes of Fig. ^ shall form the template for most of the results 
to come in the following sections. We show in Fig. f| contours of the degree of Yukawa 
unification R in per cent. The Yukawa unification varied from less than 5% to over 20% 
throughout the plane. We note that Fig. [1]g?, makes it evident that for lower values of mi/ 2 
we will have better Yukawa coupling unification occurring at slightly lower values of tan f3 
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than tan (3 = 50 in this figure. 

It was pointed out in [21| that the lightest slepton in this model is the stau, and that 
solutions with m fl < 200 GeV are hard to obtain. This is confirmed in Fig. |5| where we plot 
the stau mass as a function of the left-right mixing angle in the stau sector. Four different 
regions are indicated according to the value of the D-term mass marameter Mp (there is 
some overlap between them in the boundaries). As opposed to models with universality, 
where the light stau has mainly a right stau component, in SO(10) it can have any value for 
the mixing angle. In ISAJET, we take f\ = cos6 t ?l — sin# T f#, so that small values of 
prefer cos# r ~ and a dominantly right fi, while large values of Mp prefer cos# T ~ 1 and 
a dominantly left fx. Any intermediate value of cot# T is also allowed, so that the light stau 
could be a very mixed state of ?l and tr. 



III. EVALUATION OF NEUTRALINO RELIC DENSITY 

Recent precision measurements of cosmological parameters present new and interesting 
constraints on the possibility of particle cold dark matter p0||. Analyses of the cosmic 



microwave background radiation suggest a matter /energy density of the universe VL = = 
1.0 ± 0.2, i.e. consistent with a flat universe. In addition, the Hubble constant itself is 
now estimated as H = (65 ± 5) km/sec/Mpc, yielding a value flh 2 ~ 0.42, where H = 
lOO/i km/sec/Mpc with h the scaled Hubble constant. The baryonic contribution to Qh 2 
is estimated to be VLsh 2 ~ 0.02 from Big Bang nucleosynthesis (BBN) arguments, while 
the contribution from neutrinos ought to be Q u h 2 ^ 0.06, mainly from models of structure 
formation in the universe. Information on the total matter density of the universe can 
be found by combining constraints from BBN with measurements of intracluster gas and 
gravitational binding in rich galactic clusters: the results indicate VlMh 2 ~ 0.17. The 
remaining matter density may come from cold dark matter (CDM) particles, while the 
remaining energy density may come from a non-zero cosmological constant, as is suggested 
by recent measurements of type la supernovae. The lightest neutralino of supersymmetry 
is an excellent candidate for CDM in the universe. In a universe with 5% each of baryonic 
matter and massive neutrinos, and a cosmological constant with Q\h 2 ~ 0.25, we would 
expect Q-zJi 2 ~ 0.13. As an extreme, assuming no contribution from the cosmological 
constant, we would obtain ^^_h 2 ~ 0.38. In order to account only for the dark matter 
needed to explain galactic rotation curves, we would expect ^ 2 ~ 0.02. 

To estimate the relic density [|3T| of neutralinos in Yukawa-unified 50(10), we follow the 
calculational procedure outlined in Ref. [^]. Briefly, we evaluate all tree level neutralino 
annihilation diagrams exactly as helicity amplitudes. We then calculate the neutralino anni- 
hilation cross section, and compute the thermally averaged cross section times velocity using 
the fully relativistic formulae of Gondolo and Gelmini |33| . Once the freeze-out temperature 



is obtained via an iterative solution, we can straightforwardly obtain the neutralino relic 
density VL-^Ji 2 . Our program takes special care to integrate properly over any Breit-Wigner 

poles in s-channel annihilation diagrams p2| . We do not include Z\ — T\ co-annihilation 



diagrams; these can be important in mSUGRA models [34], but are unimportant in the 
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SO (10) model since the slepton masses are always far heavier than the lightest neutralino .0 

Our calculations of the relic density also do not include effects that might be significant |3(| 

if the 6-squark is sufficiently degenerate with Z\. 

Our results for the neutralino relic density are shown in Fig. ^ for the same parameters 

as in Fig. |. We show contours of tth 2 = 0.02, 0.1, 0.15, 0.2, 0.25, 0.3 and 0.35. In 

frame a), we see that there exists wide ranges of parameter space for which the relic density 

0.02 < ^zj 1 " 1 ~ 0.35, i.e. in the cosmologically interesting region. In no part of the 

plane shown does the relic density exceed ^^Ji 2 = 0.81. Much of the interesting region of 

parameter space occurs for large values of scalar masses, which is contrary to the situation 

at low tan f3, where large scalar masses suppress the annihilation cross section, and typically 

yield h 2 > 1. At large tan/3, the decay widths of the Higgs bosons H and A become 

very large, typically 50 — 100 GeV. Then annihilation can proceed through the very broad 

s-channel H and A resonances, even if 2m^ ^ tua or m#. In a diagonal swath that splits 

the allowed parameter space, the relic density drop to ^ z± h 2 < 0.02. This trough is due to 

regions where ttla — 2m^, and where annihilation can occur very efficiently through the 

s-channel diagrams. In frame b), we show similar results except for A = 0. Again, we see 

a region between the fl^ h 2 = 0.02 with very low relic densities. Throughout the remaining 

parameter space region, the relic density is of cosmologically interesting values, and in fact 

never exceeds h 2 ~ 0.6. 
Z\ 

IV. DIRECT DETECTION RATES FOR RELIC NEUTRALINOS 

A consequence of the SUSY dark matter hypothesis is that a non-relativistic gas of 
neutralinos fills all space. To test this hypothesis, a number of direct neutralino detection 
experiments have been built or are under construction ||37|| . The general idea behind these 
experiments is that neutralinos (or some other CDM candidate) could elastically scatter off 
nuclei in some material, transferring typically tens of keV of energy to the recoiling nucleus. 
Current generation detectors are aiming at a sensitivity of 0.1-0.01 events/kg/day 

The first step involved in a neutralino-nucleus scattering calculation is to calculate the 
effective neutralino-quark and neutralino-gluon interactions. The neutralino-quark axial 
vector interaction leads in the non-relativistic limit to a neutralino-nucleon spin-spin in- 
teraction, which involves the measured quark spin content of the nucleon. To obtain the 
neutralino-nucleus scattering cross section, a convolution with nuclear spin form factors must 
be made. The neutralino-quark and neutralino-gluon interactions (via loop diagrams) can 
also resolve into scalar and tensor components. These interactions can be converted into 
an effective scalar neutralino-nucleon interaction involving quark and gluon parton distribu- 
tion functions. A neutralino-nucleus scattering cross section can be obtained by convoluting 
with suitable scalar nuclear form factors. The final neutralino detection rate is obtained by 



2 Recently, two papers have appeared on calculating the relic density in Yukawa unified SO(W) 
p5| . Since these authors do not insist on correct third generation fermion masses, they are able 
to obtain sparticle mass spectra without invoking D-terms. Thus, their results differ significantly 
from ours. 
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multiplying by the local neutralino density (for which estimates are obtained from galaxy 
formation modeling), and appropriate functions involving the velocity distribution of relic 
neutralinos and the earth's velocity around the galactic center. 

In this section, we present event rates for direct detection of relic neutralinos left over 
from the Big Bang. For illustration, we present detailed calculations for a 73 Ge detector, 
which has a sizable nuclear spin content J = |. We follow the procedures outlined in Ref. 
P8fl , and assume a local neutralino density pj = 5 x lCr 25 gcm~ 3 [ffifl. 

Our results are shown in Fig. [7] for the same parameter plane as in Fig. |j. We show 
direct detection rate contours of 1, 0.1, 0.01 and 0.001 events/kg/day. In the low myz 
regions of parameter space, the direct detection rate exceeds 1 event/kg/day, which is a 
very large rate for neutralino dark matter. The large rate occurs generically at large tan /3 
1 4T)| , |3"8"|| . Neutralinos can scatter off quarks via squark or Higgs boson exchange graphs, and off 



gluons via loop graphs containing quarks, squarks and gluons. Neutralino-gluon scattering 
via Higgs exchange into a 6-quark loop for instance is enhanced at large tan /3, and helps lead 
to the large scattering rate. Unfortunately, the largest direct detection rates occur mainly 
in those regions of parameter space where the relic density is tiny. However, there do exist 
significant regions of parameter space where the direct neutralino detection rate exceeds 0.01 
events/kg/day where the relic density is comfortably within the cosmologically interesting 
range. 



V. CONSTRAINT FROM RADIATIVE B DECAY 

A. Results for BR(b -> si) 

The radiative flavor- changing decay b — > si can occur in the SM via Wt loop diagrams; 
in SUSY theories, additional diagrams including tH + , W{tj, gq and T^q loops will occur 



UeH- The ALEPH collaboration |43| has measured 

BR(b -> s-y) = (3.11 ± 0.80 ± 0.72) x 10" 4 (5.1) 

while the CLEO collaboration |44j has recently reported 

BR(b -> S7) = (3.15 ± 0.35 ± 0.32 ± 0.26) x 10" 4 (5.2) 

and restricts the branching ratio at 95% CL to be 

2 x 10~ 4 < BR(b -> try) < 4.5 x 10~ 4 . (5.3) 

The SM prediction at NLL accuracy j|5| is 

BR(b -> sj) = (3.28 ± 0.33) x 10" 4 . (5.4) 

The calculation of the width for b — > 57 decay proceeds by calculating the loop interaction 
for b — > si within a given model framework, e.g. the MSSM, at some high mass scale 
Q ~ Myy, and then matching to an effective theory Hamiltonian given by 

AC 8 

H eff = -Fv tb V t * s Y,d(Q)Oi(Q), (5.5) 
V2 i=1 
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where the Ci(Q) are Wilson coefficients evaluated at scale Q, and the Oi are a complete 
set of operators relevant for the process b — > 57 (given, for example, in Ref. [4q|.) All 



orders approximate QCD corrections are included via renormalization group resummation 
of leading logs (LL) which arise due to a disparity between the scale at which new physics 
enters the b — > sj loop corrections (usually taken to be Q ~ My/), and the scale at which 
the b — > S7 decay rate is evaluated (Q ~ m&). Resummation then occurs when we solve the 
renormalization group equations (RGE's) for the Wilson coefficients 

Q^C i {Q)= lji C j {Q), (5.6) 

where 7 is the 8x8 anomalous dimension matrix (ADM), and 

7 = ^ 7 (0) + (^) 2 7 (1) + • • • • (5-7) 

The matrix elements of the operators are finally calculated at a scale Q ~ m;, and multi- 
plied by the appropriately evolved Wilson coefficients to obtain the final decay amplitude. 
Recently, next-to-leading order QCD corrections have been completed for b — > 57 decay. 



These include 1) complete virtual corrections [[47| to the relevant operators O2, O7 and Og 



which, when combined with bremsstrahlung corrections |4"8]j4"?fl results in cancellation of 
associated soft and collinear singularities; ii) calculation of 0(a 2 s ) contributions to the ADM 



elements 7^ for i,j = 1 — 6 (by Ciuchini et al. |f49|), for i, j = 7, 8 by Misiak and Miinz [50 

and for 727* by Chetyrkin, Misiak and Miinz ||51|| . In addition, if two significantly different 
masses contribute to the loop amplitude, then there can already exist significant corrections 
to the Wilson coefficients at scale My/. In this case, the procedure is to create a tower of 
effective theories with which to correctly implement the RG running between the multiple 
scales involved in the problem. The relevant operator bases, Wilson coefficients and RGEs 
are given by Cho and Grinstein |)2] for the SM and by Anlauf [53] for the MSSM. The latter 



analysis includes contributions from just the tW, tH and tiWj loops (which are the most 
important ones). We include the above set of QCD improvements (with the exception of 
727^, which has been shown to be small [|51[]) into our calculations of the b — > 57 decay rate 
for the mSUGRA model. 

We include the contributions to Cj{My/) and Cg(My/) from gq and Ziq loops. To do so, 
the squark mixing matrix T which enters the couplings must be derived. To accomplish this, 
we first calculate the values of all running fermion masses in the SM at the mass scale Mz- 
From these, we derive the corresponding Yukawa couplings f u , fa and h e for each generation, 
and construct the corresponding Yukawa matrices (f u )ij, (fd)ij and (f e )ij, where i, j = 1, 2, 3 
runs over the 3 generations. We choose a basis that yields flavor diagonal matrices for (fdjij 
and {f e )ij, whereas the CKM mixing matrix creates a non-diagonal matrix (f u )ij ||27|| . The 
three Yukawa matrices are evolved within the MSSM from Q = Mz up to Q = Mqut 
and the values are stored. At Q = M GUT , the matrices (A u f u )ij, (A d f d )ij and (A e f e )ij are 
constructed (assuming A(Mqut) = A x 1). The squark and slepton mass squared matrices 
(ml)ij are also constructed, where k = Q,u,d, L and e. These matrices are assumed to be 
proportional to 1 at Q — M GUT . The (Af)ij and (ml)ij matrices are evolved along with 
the rest of the gauge/ Yukawa couplings and soft SUSY breaking terms between M z and 
Mqut iteratively via Runge-Kutta method until a stable solution is found. At Q = Mz, 
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the 6x6 G?-squark mass squared matrix is constructed. Numerical diagonalization of this 
matrix yields the squark mass mixing matrix T which is needed for computation of the gq 
and Ziq loop contributions. At this point, the Wilson coefficients Ci(M w ) and C 8 (M W ) can 
be calculated and evolved to Q ~ m b as described above, so that the b — > 37 decay rate can 



be calculated p^j55| . 



As indicated above, the main non-SM contributions to B{b — > S7) are tif* and UW. 



± and UWf. 

There is some correlation between the if ± and t\ masses and the D-term mass parameter 
M D , as we show in Fig. [8[ Large values of the charged Higgs are guaranteed if M D > 350 
GeV, which involve also large values of the top squark mass. This is the most convenient 
situation for satisfying the experimental requirements for B(b —>■ S7), because the charged 
Higgs contribution decreases. Nevertheless, in order to decrease sufficiently the chargino 
contribution, it is necessary to work in a restricted region of parameter space, as indicated 
below. In contrast to the correlation in Fig. || there is no obvious correlation between the 
charged Higgs mass and the mass of the chargino. 

Our results for the same parameter space planes as in Fig. f| are shown in Fig. |9|. 
Contours of B(b — > 37) are shown for 5, 6, 7, 8, 10 and 15 x 10~ 4 . In frame a), for 
A = — mi 6 , we see that very large values of B{b — > 57) are obtained, usually well above the 
CLEO 95% CL limit, except for the very largest values of mi/ 2 ~ 1000 GeV. In particular, 
the regions with the maximum rate for direct detection of dark matter also have B(b — ► 57) 
values far in excess of experimental measurements. A similar behavior is seen in frame b) for 
A = 0. The very large values of B{b — > 57) are typical of SUSY models with large values of 
tan/? and [i < 0. The qualitative behavior that Yukawa unification works best for the sign 
of fj, which gives the largest values for B(b —> 57) has been noted previously by a number of 



authors 14 . 5C 



We do note here that there may be other possible effects which could act to reduce 
the b — > 57 branching fraction. Nontrivial flavor mixing in the down squark sector could 
result in significant non-cancelling contributions from gluino loop diagrams, thus reducing 
the value of B(b — > 57). Furthermore, additional sources of non-universality could act to 
yield values of B(b — > sj) more nearly in line with experimental values ||56|| . For instance, 
non-universality between generations is still allowed by 50(10), and could lead to potentially 
large gluino loop contributions. This could, for instance, naturally be the case in a class 
of 5*0(10) models where dynamics makes third generation sfermion masses much smaller 
than masses of other sfermions |57|]. Also, modifications can occur by including large CP 



violating phases for terms such as the /1 parameter which controls the alignment of the SM 
Wilson coefficients with the top squark contribution. Of course, such phases would also 
affect for instance the electron and neutron electric dipole moments, and even scattering 



cross sections |p8| . Finally, radiative corrections due to Yukawa couplings are not included 
in the computation of the decay rate. Keeping these issues in mind, one should exercise 
judgement in drawing broad conclusions from the results of Fig. |]. 

It is well known that the supersymmetric contribution to B(b — > 57) is proportional to 
A t /j,ta,nf3 when tan/3 is large |5!J. Motivated by this, we have made a dedicated search in 
a restricted area of parameter space of our specific model where A t is close to zero at the 
weak scale. To be more specific, we consider: 

1300 GeV < mie < 1700 GeV 
1300 GeV < m w < 1700 GeV 
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300 GeV < M D < 600 GeV (5.8) 

50 GeV < mi/a < 300 GeV 
800 GeV < A < 1800 GeV 

In Fig. [ED), we plot the branching ratio B(b —>■ 37) as a function of the light chargino mass, 
including only the points that satisfy the experimental constraints from colliders (aside from 
the LEP2 chargino mass bound). We see that many solutions can occur with B(b —>■ 37) 
within the 95% CL range, and with a low value of m^. These small A t solutions typically 
have 777.10 < m m, and very low chargino masses, but top squark masses above 1 TeV, which 
acts to suppress UWi loop contributions. They correspond to solutions mostly to the left of 
the diagonal line in the upper part of Fig. [j]a). 

B. Comparison with SO(10) analysis of Blazek and Raby 

A recent study by Blazek and Raby has appeared |JD| (BR), which also focusses on the 
rate for b — > 37 decays in Yukawa unified 50(10) models. There are a number of differences 
between their study and ours. The main difference is that BR, based on a global fit [^TJ of 



model 4c of Anderson et al. (Ref Q), reject all models with fi < based on the generally 
large rate for BR(b — > 57) obtained. On the contrary, we reject all models with fi > 
based on the inability to achieve a high degree of Yukawa coupling unification. BR adopt a 
top-down solution to renormalization group running, and impose universality of all matter 
scalars at Q — M GUT , although they allow for arbitrary GUT scale values of m Hu and mn d , 
whose values are fixed by imposing REWSB. This allows BR to generate solutions with very 
small values of /1, so that their Z\ will be largely Higgsino-like. They also generate rather 
small SUSY loop corrections to mj, ranging from 0% — 8%. These may be in part due to 
the small values of a s [Mz) they generate, ranging from 0.11 — 0.18. In contrast, for fi > 0, 
we generate SUSY loop corrections to m b of order 10% — 40% at large tan j3, in accord with 
Pierce et al., Ref. ||22|1 . Using a bottom- up approach to renormalization group running, the 
large SUSY loop corrections to fermion masses allow us to generate Yukawa unification only 
good to 25% — 75% for this sign of /i (again in accord with Pierce et al). Thus, if we allow 
a much more liberal constraint on Yukawa unification, we will also obtain fi > solutions, 
which in general can give BR(b — ► 37) rates more closely in accord with measurements. We 
note that for /j, > 0, we find REWSB can occur for tan/5 ~ 50, even assuming universality 
of all scalar masses (as mentioned earlier in Sec. I). 



C. Note added: 



After completion of this manuscript, a preprint appeared p2| wherein it was claimed 



that NLO corrections to the b — > sj loop diagrams can lead to sign reversals on 

certain amplitudes, especially at large tan/5. This effect leads to rates for b — > 57 decay 
in accord with theory predictions for Yukawa unified models with \i < 0, while excluding 
models at large tan (3 with \i > 0. In view of these (incomplete) calculations, it appears to us 
that one should view the restrictions from one-loop calculations of b — > sj decay with care. 
Rather than categorically excluding regions, it appears to be best to wait until complete 
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calculations are available. It would indeed be exciting if these corrections led to values of 
the 6 — > 57 decay rate in agreement with experiment for the sign of mu favored by Yukawa 
unification. 



VI. PROSPECTS FOR YUKAWA-UNIFIED 50(10) AT THE FERMILAB 

TEVATRON 

In Yukawa unified SO (10) models, as with most SUSY models, the light Higgs boson is 
bounded by < 125 GeV, and so may well be accessible to Tevatron searches with high 



luminosity and improved detectors |65|j66|j59|1 . The h boson will in general be difficult to 



distinguish from a SM Higgs boson, so that signatures from other sparticles or additional 
Higgs bosons would be useful. In much of the parameter space of Yukawa unified 5*0(10) 
models, the gluinos, squarks and sleptons are too heavy to result in observable signals for 
the Fermilab Tevatron collider. There are, however, three exceptions to this. 1.) For low 
values of mig, the pseudoscalar and heavy scalar Higgs bosons A and H may be accessible 
to searches via pp — > bbA or bbH. 2.) If tan/3 is on the lower end of acceptable values 
(e.g. tan/3 ~ 46 — 49), then some regions of parameter space with low values of mu2 may 
be accessible to trilepton searches. 3.) For large negative values of the A parameter, the 
&i-squark becomes light enough that there exists a substantial production cross section for 
&i&i production. In this section, we address each of these possibilities.^ 



A. Search for the H and A Higgs bosons 

Here, the search is for bb(f> events, where <fr — A or H. At large tan /3 ~ 50 as in Yukawa 
unified S'O(IO), the bbA and bbH couplings are enhanced, leading to large production cross 
sections. For the Higgs mass range of interest, H and A typically decay to 66 with a ~ 90% 
branching fraction, so that the signal will be four 6-jet events, where two of the 6s reconstruct 
to the Higgs mass. For Higgs masses in the range of interest, usually m# ~ ttla, so that the 
resonances from each will be overlapping. The dominant backgrounds come from 6666, bbjj, 
Wbb and Zbb production. 

In Fig. [11] we plot allowed regions in the m w -mA plane when the independent parameters 
of the theory are varied as indicated in Eq. fl2.ip . Four regions are displayed according to 
the value of the D-term mass parameter Mo (the four regions have a small overlap in their 
boundaries). The following correlation is obvious from the figure: a smaller CP-odd Higgs 
mass rriA is obtained for smaller and mw mass parameters. We have checked that the 
heaviest CP-even Higgs mass m# is very close to rriA, in fact they are within ±2% (±5%) 
of each other for > 200 (100) GeV. On the contrary, if ijia is decreased below 100 GeV, 
H becomes heavier than A, such that mu can be up to 50% larger than itia for itla ~ 80. 



3 The effect of SO (10) D-terms on SUSY cross sections at the Fermilab Tevatron collider has been 
examined in Ref. j67|. These authors do not insist on Yukawa unification, so that any value of 
tan p is allowed. 
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For the recently concluded year-long Tevatron workshop on SUSY/Higgs physics at Run 
2, analyses have been performed to estimate the reach for H and A at large tan/3 for both 
the DO and CDF detectors [P^j . In this report, details of expected event rates, backgrounds 
and selection cuts can be found. For tan (3 = 50, the 5a discovery reach of the DO detector 
is expected to be uia — 145, 190 and 220 GeV for integrated luminosities of 2, 10 and 30 
fb" 1 , respectively. Similarly, the reach of the CDF detector is estimated to be rriA — 130, 
160 and 220 GeV for the same integrated luminosities. 

In Fig. [12], we show contours of tua that correspond to the 5a reach of the DO detector, 
for the same parameter choices as in Fig. |], assuming integrated luminosities of 2, 10 and 
30 fb -1 . The region to the left of the contours ought to be accessible to Tevatron searches 
for bbH and bbA events. These regions represent the most likely possibility for detection of 
beyond the SM physics at the Tevatron if Yukawa unified S'O(IO) is indeed correct. 



B. Search for trilepton events 

One of the most promising ways to detect supersymmetry at the Tevatron for models 
with universality and low values of tan/3 is via trilepton events, which dominantly come 



from W\Z 2 — > 3£+ $t p8| . At large values of tan/3 and low values for scalar masses, the 
branching fraction of Z 2 is dominated by real or virtual staus, so that the Z 2 branching 
fraction to es or /is is suppressed. Models with < 200 GeV have the potential to yield 
significant trilepton signals at the Fermilab Tevatron pp collider. 

In a recent paper |69|] , the trilepton signal at the Fermilab Tevatron was re-examined for 



the mSUGRA model, in light of new backgrounds from W*Z* and W*"f* production and 
decay, and the inclusion of decay matrix elements into ISAJET. New sets of soft (SC2) and 
hard (HC2) cuts were proposed which allowed the trilepton signal to be seen at the Tevatron 
collider for integrated luminosities of 2 and 25 fb -1 . For the Yukawa unified 5*0(10) model, 
we have used the background rates and the SC2 cuts presented in Ref. [^] to evaluate the 
observability of various points in model parameter space. 

In scans of the parameter space corresponding to Fig. |, no points were found where 
the Tevatron could see the 3£ signal using cuts SC2 with either 2 or 25 fb -1 of integrated 
luminosity. For both frames a) and b), the bottom squark b\ is relatively light compared to 
other scalars so that Z 2 dominantly decays to bbZi final states, at the expense of e and /1 
modes. 

By modifying the parameter space somewhat, it is possible to find significant regions 



where the 3£ signal should be visible. In Fig. [I3, we show the m 16 vs. mx/2 plane for 
m 10 = 1.14m 16 , A Q = 0, tan/3 = 47 and = ^m^. The mass contours shown are the 
same as in Fig. ^|. In this case, there is a significant region below the =150 GeV 
contour accessible to Tevatron searches. No points were found to be observable at 5a with 
just 2 fb -1 of integrated luminosity, using cuts SC2. However, the open squares denote 
points where a 5a 3£ signal can be seen above background with 25 fb -1 . The diamonds 
denote points seeable at 3a level with 25 fb -1 of data. In the region of parameter space 
with Tn,yp < 150 GeV but with low values of m^, the 61 is becoming sufficiently light that 

three body Z 2 decays to 6-quarks dominates. In the visible region at higher values of mi6, 
the squarks, including b\ are so heavy that the Z 2 decay is dominated by the virtual Z 
contribution, and Z 2 — > e + e~Z\ with typically a 3% branching fraction. If we use the same 
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parameters as in Fig. [13] but increase tan f3 to 50, then the b\ becomes so light that again 
Z 2 dominantly decays to 6-quarks, and the trilepton signal vanishes. 



C. Search for direct production of 6-squarks 



In the Yukawa unified SO (10) model, it has been noted in Ref. |2T] that a combination 
of the large 6-quark Yukawa coupling plus the .D-term contribution to the b^ squared mass 
yields parameter ranges with the light 6-squark b\ being by far the lightest of all the squarks. 
In fact, frequently b\ is also much lighter than all the sleptons. This means that there may 
be regions of model parameter space where b\b\ pairs can be produced at large rates, with 
possibly observable signals. 

In Fig. H we show the lightest bottom squark mass m~ b as a function of mi6 for different 
values of the D-term parameter Mb- The three regions correspond to Mb < 220 GeV, 
220 GeV < Mb < 330 GeV, and 330 < Mb GeV, with some overlap at the boundaries. The 
smallest bottom squark masses we obtain satisfy < 200 GeV and occur for intermediate 
values of Mb and m^. These points have approximately /x ~ —200 GeV imposed by EWSB 
(and Aq ~ —1600 GeV). This value of \x together with the large value of tan/3, makes the 
left-right sbottom mixing comparable with the right sbottom soft mass parameter mo. 

In Ref. |[29|| , the reach of the Tevatron collider for b\b\ pairs was investigated. If b\ — > bZ\ 
is the dominant decay mode, then the final state will include two hard fe-jets plus fix- 
Backgrounds investigated included W+jets, Z+jets and tt production. For the vector boson 
production backgrounds, the jets qualifier includes both 6-jets and other light quark jets. A 
number of cuts were proposed, including tagging the 6-jets via microvertex detectors, which 
allowed signals to be seen above backgrounds in large parts of the vs. parameter 
space. Typically, bottom squark masses of m~ bi ~ 210 GeV could be seen with 2 fb -1 , and 
mi ~ 240 GeV could be seen with 25 fb -1 . The reach was somewhat reduced if b\ — > bZ 2 
decays were also linematically accessible. 

In Fig. [1^, we show the same parameter space as in Fig. |]a). We have generated 
events with ISAJET for pp collisions at y/s = 2 TeV using the detector simulation, cuts and 
background rates described in Ref. PP[. Parameter space points labelled with a star are 



seeable at the 5a level with just 2 fb 1 of integrated luminosity. In addition, points labelled 
with an open square are seeable at 25 fb -1 . The open squares have b\ masses as heavy as 
245 GeV. 



VII. SUMMARY, CONCLUSIONS AND OUTLOOK 

We have shown before that supersymmetric particle spectra can be calculated in minimal 
super symmetric 50(10) GUT models consistent with Yukawa coupling unification to 5% or 
better, a top mass of mt = 175 GeV and including radiative electroweak symmetry breaking. 
To accomplish this, we made crucial use of D-term contributions to scalar masses that are 
induced by the reduction in rank of gauge symmetry when 5O(10) breaks to the MSSM at 
Q = Mqut- The D-terms split the GUT scale values of the Higgs boson squared masses so 
that m% u < rn 2 Hd already at Mqut- This allows Tn 2 Hu to be driven more negative than rn 2 Hd 
at Q = M wea k, which is required for REWSB to take place. We included the essential SUSY 
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loop corrections to third generation fermion masses to generate the appropriate weak scale 
Yukawa couplings. Yukawa unification to 5% takes place only for /j, < and tan/3 ~ 46 — 52. 
In this paper, we upgraded our Yukawa coupling evolution equations to include two-loop 
terms in the RGEs plus weak-scale threshold effects for the one-loop terms. 

The sparticle mass spectra we calculate reflect the influence of the Z)-terms, and as 
a consequence, the b\ turns out to be almost always the lightest of the squarks, and is 
sometimes even lighter than the sleptons. In addition, the lightest sleptons can frequently 
be (dominantly) left- sparticles, as opposed to the mSUGRA model where they are usually 
right- sparticles. This latter trait may be tested explicitly at e + e~ linear colliders operating 
above slepton pair production threshold. 

We attempted to systematize the parameter space of these models, showing plots of 
sparticle mass in the m\§ vs. m 1( / 2 plane for viable selections of mio, Md and tan/3. Of 
course, fixing these parameters meant we had to loosen our restrictions on the degree of 
Yukawa coupling unification. We stress though that the original degree of unification can 
be regained by a small change in tan /3. 

We calculated the relic density of neutralinos in Yukawa unified 5*0(10), and found none 
of the parameter space that was explored to be excluded by the limits on the age of the 
universe. In fact, very desirable values for Q^Ji 2 were found throughout much of parameter 
space. 

We also evaluated the direct neutralino detection rate for a 73 Ge detector. Detection rates 
could exceed 1 event/kg/day, but only in parameter regions with very low relic densities. In 
more favorable regions, the direct detection rate often exceeded 0.01 events/kg/day, which 
is the capability being sought by current generation dark matter detection experiments. 

We also examined the branching ratio for b — > 57 decay in Yukawa unified 50(10). 
We found that it was beyond the 95%CL experimental limits for much of the parameter 
space shown. However, special regions of parameter space giving rise to small values of 
weak scale A t gave more acceptable B(b — > sj) values, but with > 1 TeV. There 
could be other mechanisms such as non-trivial mixing in the squark sector, additional non- 
universality or mixing, or CP violating phases could act to reduce the large B(b —>■ 37) 
values. Alternatively, the rate for b — > 37 is more favorable for \i > 0, but then the degree 
of Yukawa unification would have to be relaxed substantially. Finally, recent claims have 
been made |^2j that two-loop b — > S7 calculations can actually reverse the signs of certain 
decay amplitudes, leading to the possibility of agreement between b — > 37 calculations and 
Yukawa unified SUSY models with \x < 0. In light of these claims, it seems prudent to 
refrain from excluding parameter space regions based on one-loop b — > 57 calculations until 
the theoretical predictions become more certain. 

We found three possibilities for new particle detection at the upgraded Fermilab Tevatron 
beyond discovery of the lightest Higgs boson. In regions of parameter space with small m\e, 
bbA and bbH production could be visible. In some parameter space regions with small m^, 
it is possible to detect trilepton signals from W\Zi — > 3£ production. In other regions again 
with small mi/2, bottom squark pair production may yield a visible signal. The other scalar 
particles are generally too heavy to be produced substantially at the Tevatron collider. 

At the LHC, we expect Yukawa unified 5O(10) to be visible as usual for values of rrig 
beyond ~ 2 TeV. Within this framework, SUSY events ought to be rich in 6-quark jets due 
to the relatively light third generation squark masses inherent in this theory. 
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Yukawa unified SO (10) models have been historically very compelling models for reasons 
listed in the introduction, and are even more compelling today due to recent strong evidence 
on neutrino masses. The Yukawa unified SO (10) model is only one of a number of SO(10) 
models that are possible. A more general SO (10) model p4] could contain a superpotential 
with terms such as 

f = hUid + ftUk + ''' (7-1) 

in which case just b — r unification would occur. This would allow a much wider range of 



tan p ~ 30 - 50 to occur |24] 



In addition, recent work on radiatively driven inverted hierarchy models has shown that 
"natural" models can be obtained with multi-TeV scale scalar masses, where the third 



generation masses are driven to weak scale values |70]| . These models suffer the same problem 
with REWSB (and possible color and charge breaking minima in the scalar potential) as 
the SO(10) model considered here. It has recently been shown that D-term contributions to 



scalar masses aid this class of models to acheive REWSB | 57| . Finally, more complex gauge 
symmetry breaking patterns could occur which would affect scalar masses in alternative 
ways 0. 
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FIG. 1. Plots of regions of parameter space where valid solutions to minimal SUSY SO(10) are 
obtained, consistent with Yukawa coupling unification to 5%, and radiative electroweak symmetry 
breaking. We use 2-loop RGEs for both gauge and Yukawa coupling evolution. 
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FIG. 2. We plot theoretically allowed regions of m\% vs. parameter space taking Aq = 0, 
mio = 1-25tt7i6, tan/3 = 50 and \i < 0. In frame a) we take Mp = \m\§, while frame b we take 
Md = |tt7i6. The shaded areas are not allowed by the REWSB constraint. Open dots represent 
experimentally excluded points. The dashed contours, from bottom to top, are for m„ L = 1000 
and 2000 GeV. The solid contours, from bottom to top, represent masses of 150, 250 and 

350 GeV. The dotted contours represent from left to right txia = 150, 250 and 350 GeV. Yukawa 
coupling unification varies from less than 5% in the central regions of the plot, to nearly 25% at 
the edges of the excluded region. 
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FIG. 3. The same as Fig. ||, except Mc = |mi6, while in a), Aq = —miQ and in b), Ao = +mie. 
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FIG. 4. Contours showing Yukawa coupling unification in the m\% vs. m 1 / 2 plane, for 
Mr> = |mi6, tan/? = 50, mio = 1.25mi6, A* < and a) ^4o = —mw and 6) = 0. In frame b) the 
gree of unification reduces to 5-10% beyond the right-most contour in the figure. 



26 




0.0 0.2 0.4 0.6 0.8 1.0 



COS 0y 

FIG. 5. Lightest stau mass as a function of the stau mixing parameter cos 2 9 T for different 
values of the D-term parameter Mr>. In 5O(10) the stau mixing can take any value. 
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FIG. 6. The same parameter space planes as in Fig. £§ but now showing contours of neutralino 
relic density % h 2 = 0.02, 0.1, 0.15, 0.2, 0.25, 0.30 and 0.35. 
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a) A..=-m., b) A ft =0.0 
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FIG. 7. The same parameter space planes as in Fig. |], but now showing contours of direct 
neutralino detection rate in events/kg/day in a 73 Ge detector. 
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FIG. 8. Charged Higgs mass as a function of the light top squark mass for different values of 
the D-term parameter Mp. Heavy and t± are associated with large values of Mp, which is 
favoured by B(b ->■ S7). Here, /x < and \A \ < 3 TeV. 
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a) A..=-m., b) A ft =0.0 
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FIG. 9. The same parameter space planes as in Fig. |4|, but now showing contours of decay 
branching fraction for b — ► 57. Starting from the top, the contours are for a branching fraction of 
5,6,7,8,10 and 15 times 10" 4 . 
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FIG. 10. Branching ratio B(b — > sj) as a function of the light chargino mass when the input 
parameters are varied within the intervals given in Eq. ( |2.1| ). These points in parameter space are 
characterized by small \At\, which favours small values of B(b — ► s'y) when tan/3 is large. 
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FIG. 11. Neutral CP-odd Higgs mass vtla as a function of mio for different values of the 
D-term parameter Mp. The heavy neutral CP-even Higgs H has a mass very close to rriA except 
when rriA < 100 GeV. 
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FIG. 12. The same parameter space planes as in Fig. ||, but now showing regions of parameter 
plane where 5a bbA and bbH signals may be seen with the DO detector at Fermilab Tevatron for 
(contours from left to right) integrated luminosities of 2, 10 and 30 fb . 
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FIG. 13. Regions of the m\% vs. parameter plane where trilepton signals may be seen 

at the Fermilab Tevatron. Open squares can be seen at the 5a level with 25 fb _1 of data, while 
diamonds can be seen at the 3a level. We take Aq = 0, tan/3 = 47, mio = 1.14mi6, Mu = 
and fi < 0. 
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FIG. 14. Light bottom squark mass as a function of m\% for different values of the D-term 
parameter Md- Very light sbottoms can be obtained for intermediate values of Md- 
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FIG. 15. Regions of the mi6 vs. m 1 / 2 parameter plane where b\bi — > 66+ events can be 
seen above background at the Fermilab Tevatron. The stars show points where a 5a signal can be 
seen at Run 2 (2 fb _1 ) and the open squares show points where a 5a signal can be seen at Run 3 
(25 fb _1 ). We take Aq = —mie, tan/3 = 50, m\o = 1.25mi6, = |mi6 and /i < 0. 
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